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We calculate the tunneling current in a quantum Hall bilayer system in the strong mag- 
netic field limit. We model the bilayer electron system as two Wigner crystals coupled 
through interlayer Coulomb interactions, treated in the continuum limit. We generalized 
the Johansson and Kinaret (JK) model^ and were able to study the effect of the low 
energy out-of-phase magnetophonon modes produced as a result of tunneling events. We 
find the same scaling behavior of the tunneling current peak with the magnetic field as 
found by JK but were able to find the tunneling current scaling behavior with interlayer 
distance as well. 
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1. Introduction 

A two-dimensional bilayer electron system under the presence of a strong perpen- 
dicular magnetic field exhibits a Coulomb barrier peak in the tunneling current in 
transport experiments^. This peak is pronounced in high density samples and is 
attributed to the highly correlated state of the electron liquid where the cyclotron 
motion frustrates the ability of local electrons to respond to a tunneling event. A 
necessary amount of energy of the order of e^/47re(a), where (a) is the average in- 
terelectron distance, has to be provided to the system to collectively respond to 
tunneling. 

Our model describes this tunneling behavior and is based on the Johansson 
and Kinaret (JK) model^ which assumes that the two layers of electrons are in 
commensurate Wigner crystal states. A tunneling event disrupts the Wigner crystals 
producing collective fluctuations (magnetophonons) that dissipate the energy of the 
tunneling electron and smooth-out the local charge-defect associated with it. We 
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work in the continuum limit treating the whole system as an elastic medium. We 
should stress out that the actual state of the system is not a Wigner crystal but 
a rather highly correlated liquid state. Our work differs from the JK model in the 
fact that we have added interlayer Coulomb interaction, correlating the two layers, 
and additionally we were able to produce analytic results as well as numerical ones. 
Our results are summarized in Fig. (1). 



2. Bilayer theory 



We assume that each layer of electrons is a two-dimensional elastic medium in- 
teracting via Coulomb interaction (intralayer). Each electron has a displacement 
associated with its equilibrium position in the Wigner crystal which becomes the 
displacement field in the continuum limit. The magnetic field couples the transverse 
and longitudinal modes into gapless magnetophonon and gapped (by the cyclotron 
frequency) magnetoplasmon modes. These correspond to the internal modes of each 
layer when they are far apart ( JK model) . When the layers arc brought within prox- 
imity the interlayer Coulomb interaction will introduce an additional coupling of 
those four modes producing two in-phase and two out-of-phase modes. The in-phase 
modes resemble single layer modes but the out-of-phase modes are gapped by the 
energy cost to shear the commensurate Wigner crystals. This energy cost is asso- 
ciated with the penalty the short range part of the interlayer Coulomb interaction 
introduces to an out-of-phase motion of the two Wigner crystals. Additionally, a 
tunneling electron couples only to out-of-phase mod(^s. This happens because during 
tunneling the positive charge defect the electron leaves behind has to be "covered- 
up" while the receiving layer has to create a place for the incoming electron by 
"opening-up" resulting in a relative out-of-phase motion of the two crystals. The 
total Hamiltonian for the bilayer system assumes the form, 

H = Hq -\- Hrp + 



eA + i^MsA{al 



as) 



c\cA + 



eB+i^MsB{al 



cLcb 



(1) 



where e^, cb are the corresponding Madelung energies for the two Wigner crystals, 
al, as are the creation and annihilation operators for the collective electron modes 
while c^, c are the tunneling electron's creation and annihilation operators for the 
two layers. The ground state of this Hamiltonian is a bath of collective modes 
coupled to an external electron present in the system. Tunneling is an independent 
single particle process where the tunneling matrix element T is assumed to be real 
and calculated as in JK^. The tunneling current for different bias voltage V is given 
by 
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dte-^*/ft([i/^(i),i7+(0)]), 



(2) 
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We calculate the zero-temperature limit for the correlation function, since the ex- 
perimental temperature range is O.IK 10~^ eV while the applied bias voltage is in 

the mV range. Additionally wc concentrate on the magnetophonon modes only that 
are responsible for the low bias Coulomb barrier peak. We find that the correlation 
function (iJ^ (t)i?^(0)) ^ C{t) can be Fourier transformed and obeys an integral 
equation of the form 

u>C{uj)^ [ dxf{x)uj{x)C{Lj~Lu{x)), (3) 



where uj{x) is the magnetophonon dispersion relation and f{x) is a weight function. 
Equation (3) has an approximate large frequency solution of the form 

CM^expl-l^i^l, (4) 



2C2 

where ^/c^ I/oq^/B for layers far apart and C2 ^ d/aQ^/B for layers close together, 
where d, is the interlaycr distance and hq the Wigner crystal lattice parameter. Based 
on that we propose a trial solution of the form 

C{lo) = Nu'^e-^^^ (5) 
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where N, r, A are self-consistently evaluated from the properties of C{t) namely, 

o 

dujC{uj) = 2n (6) 

duicoC {(jj) = 27rci (7) 

dcjuj'^Ciuj) = 27r(c2 + cj). (8) 
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Our solutions are presented in Fig. (1) for different applied magnetic field values. 

For all our model parameters we have used the experimental system reported in 
Ref. 2. For the magnetophonon gap value we extracted it from time-dependent 
Hartree Fock calculations of the magnetophonon dispersion relation for the bilayer 
system^. We numerically integrated Eq. (3) to verify the validity of our analytic 
solution with very satisfactory results. 



3. Conclusions 

We calculated the tunneling current curve for a bilayer quantum Hall system under 
the presence of a strong perpendicular magnetic field and reproduced the Coulomb 
barrier peak experimentally found. Our model assumed that the two dimensional 
electron systems were in Wigner crystal states and calciilatcd the eigenmodes in 
the continuum approximation where the two layers are treated as elastic media. We 
found that by introducing an interlayer Coulomb interaction the internal modes of 
the system couple to in-phase and out-of-phase ones. The latter acquire a gap in the 
magnetophonon branch associated with the short-range interlayer correlations built 
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into the system by the Coulomb force. This magnetophonon branch provides the 
necessary coUective mode couphng for tunnchng electrons to dissipate their energy, 
and for the system to relax the charge defect a tunneling event creates. We found 
the same scaling behavior JK report for the tunneling curve but additionally were 
able to extract the dependence on interlayer distance as well. 
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